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Abstract 

We investigate the transport properties of a model of an interacting electronic res- 
onant level system. The hybridization of the localized resonant level is treated in 
perturbation and the contributions to all orders are computed. This includes an 
exact treatment of the electron scattering from the localized level. The renormal- 
ization of the direct tunneling between the electrodes arising from the Coulomb 
repulsion of trapped particles is taken into account. Further, we provide a detailed 
solution of the system of coupled singular integral equations which determine the 
propagator in the interacting system. 



1 Introduction 



Over 30 years ago, Nozieres and De Dominicis (ND) [|^] proposed an ingenious 
solution of the x-ray absorption/emission edge singularity problem in metals. 
Their solution, being non-perturbative, allowed to circumvent the difficulties 
arising from the Anderson orthogonality catastrophe and revealed the sin- 
gular behavior of the absorption/emission resonance. Soon thereafter, Yuval 
and Anderson used the formal resemblance between the x-ray problem and 
the Kondo spin 1/2 model to extend the method to describe the vacuum-to- 
vacuum amplitude (or the partition function in the finite temperature case). 
They showed that the latter is formally equivalent to the partition function 
of a classical one-dimensional gas of particles with logarithmic interactions. 
In a subsequent work, Anderson, Yuval and Hamann 0] introduced renor- 
malization group methods to solve this partition function. Another formal 
correspondence to the x-ray problem can be found in the description of the 
tunneling of an electron through a localized state in a junction. For such a sys- 
tem, Matveev and Larkin |^ applied ND's technique to express the resonant 
behavior of the tunneling current. 
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The common property shared among these examples is a time- dependent in- 
teraction between conduction electrons in certain quantum channels with a 
localized structureless scatterer (a hole, an impurity spin, or a trapped elec- 
tron). Since the interaction does not change the internal state of the scatterer, 
a reduction of the many-body problem into one-body problems is possible. 
The latter can then be solved by the methods devised by ND. 

In the examples above it was assumed that electrons in different quantum 
channels (such as angular momentum, spin, pseudo spin, or indexing spatially 
separated electrodes) do not interact. Generally this results from symmetry 
arguments and a partial diagonalization of the Hamiltonian. In the present 
paper, however, we focus on the case in which such a decoupling is not possi- 
ble, and in which particles can be scattered onto different channels due to the 
interaction with the localized scatterer. We choose as an example the physical 
system considered by Matveev and Larkin of two electrodes connected to a 
barrier containing one localized state. However, in contrast to Matveev and 
Larkin, we assume a direct coupling of the conduction electrons in both elec- 
trodes. The method is illustrated through the computation of the tunneling 
current, expressed in the form of an infinite series (see Eqs. (59) and (60) 
which are the results of this paper), although a large part of the derivations 
is performed in the most general case accessible by this approach. 

This is a technical paper and we focus mainly on the formal aspects of the 
problem. We expect to come back in the future to physical applications of 
the generic problem solved below. On the other hand, the physical nature of 
the system considered below has already been discussed thoroughly in . For 
pedagogical reasons details of the computation are written down explicitly. In 
particular, we present in the Appendix A a full solution of the singular integral 
equations determining the one-particle propagators. In this respect, we also 
restrict ourselves to an equilibrium description of the system. An extension to 
the non-equilibrium case will be performed in a subsequent work. 

In Sec. 2 and 3 we define the model and derive the Dyson series for the 
tunneling current. Sec. 4 and 5 contain the central part of the paper in which 
first ND's method is adapted to the present problem, and then extended to 
all orders in the Dyson series following Yuval's and Anderson's ideas. The 
application to some simpler problems is discussed in Sec. 6. 



2 Model 

The Hamiltonian of the system reads 

H = Ho + Hw + Hv, (1) 



2 



with 



Ho = J2^kaCl^Cka + EddU, (2) 
ka 

Hw=Y: W^,i ClCk>a'd^d, (3) 

kk"a' 

Hv^^{v,"Ck^d^ + h.c.), (4) 

ka 

where and Cka are the creation and annihilation operators of spinless 
fermions with momentum k which are characterized by an additional quantum 
number a — 1,. . . ,N for some integer A^. These particles are assumed to 
form for every a a Fermi liquid with chemical potential fia, described by Hq. 
The Hamiltonian Hw describes the interaction of the particles {k, a) with a 
localized state, 0^, which is characterized by the operators d^ , d and the energy 
Ed. The hermiticity of Hw imposes that {W^,^')* = W^,'^. Finally, Hy denotes 
the hybridization term with the localized state. 

It is essential for the following to assume that this localized state is completely 
structureless and that the interaction with the Fermi liquid is short-ranged. 
Under these assumptions the short-ranged potentials can be considered to be 
separable as 

w^,i = w--'{utrut:, (5) 

Vr = \/"<, (6) 

with u'^ some normalized functions that arc peaked at the chemical potentials 
/Xq, where they take the value = 1. It is convenient to define the operators 

Ca^T^Ckaut and Ct=5^cLKr (7) 

k k 

that annihilate and create a particle localized near (p^. The Hamiltonians H^r 
and Hv become 

Hw^J2W'"''ClC„,dU, (8) 

aa' 

Hv^Y.{V"Cad^ + h.c.). (9) 

a 

As an example for such a system we consider the situation depicted in Fig. 1. 
An insulating barrier is assumed to be connected on two sides to electrodes 
containing spinless Fermi hquids. The index a labels here the left or right 
electrode, a — L,R. The localized state (pd corresponds to an impurity trap 
inside the barrier. Both electrodes are considered as particle reservoirs that 
are kept at the chemical potentials fiL and fi^ (with generally fiL 7^ /U^), 
respectively. It is assumed that in absence of a particle in the impurity trap the 
two Fermi gases are completely independent (we neglect the direct tunneling 
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Fig. 1. Sketch of the physical system underlying the model. 

which can be regarded as a constant background current). However, if an 
electron is trapped at the impurity site, it interacts with the particles near 
the junction in the electrodes since the potential created by its charge is not 
screened in the insulator. We denote this scattering potential by Wj^ff where 
the diagonal terms with respect to a, a = a', describe the backscattering 
of electrons and the non-diagonal terms, a ^ a', the renormalization of the 
background tunneling current. 

We intend to calculate the current through the junction, which can be ex- 
pressed as the time variation of the number of particles in the electrodes. 



J{t)^e(^^^{N^it)-N\t))y 



(10) 



where N°'{t) — J2kClai't)^kait) is the electron number operator in the elec- 
trode a. Its time derivative is determined by the Heisenberg equation 



(11) 



where we have defined 



N^{t) = {-t)[N''{t),Hw{t) 

a' 

N^{t) = {-t)[N-{t),Hv{t)_ 

= H){v-cimt)-ivyc^{t)S{t)}. 

Therefore the computation of the averages 



[12) 



(13) 



{N^it)) = 23m 



j:w-' {ciit)cAt)dHt)dit)) 



{m))^2jm[v^{cimt)) 



(14) 
(15) 
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which can be expressed by the response functions, 



F^'^'it) = {Cl{t)CAt)dHt)d{t)) 



(16) 
(17) 



is required. 



3 Dyson series 



We expand the functions 



F^' (t) and Fy(t) in a zero-temperature perturbation 



series in powers of Hy. The interactions due to Hw will be treated in the exact 
way indicated by ND. For this purpose, we define an interaction picture with 
respect to Hi = Hq + Hy/ by 



for an arbitrary Schrodinger operator A. In order to take averages over op- 
erators we impose artificially to the system to be in the same state \i) at an 
initial time to, as well as at the final time ti. We choose \i) to be an eigenstate 
of Hi in which the state ipd is unoccupied, d\i) = 0. This choice of the aver- 
age is definitely wrong because the occupancy of 0^ in the final state is not 
well-defined due to the hybridization term Hy. A more subtle investigation 
of the problem would involve Keldysh's formulation of perturbation theory 
for non-equilibrium processes P] in which only the initial, but no final state 
appears. An extension to the Keldysh formulation of the methods described 
below will be the subject of a future publication. In the present paper, how- 
ever, we choose for pedagogical reasons to restrict ourselves to the case in 
which the system comes back to the state \i) at the time ti. 

The standard perturbation theory tells us that in this case the average value 
over the operator A{t), for < t < ti, can be expressed by the relation 



1 



(18) 



{A{t)) = {t\T{S{to,ti)A{t)}\z) 



where T denotes the time-order operator, and where 




(19) 



is the ^'-matrix. In this representation, the functions F{^' (t) and Fy{t) read 



F^^'it) = ^\T[s{toM)&MCo.'{t)d^{t)d{t)}\i), 
F^{t) = {i\T[S{toM)Ci{t)d{t)}\i), 
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or, explicitly, 

Fw{t) = T.7^ dn...dr2„ (20) 
(z|T{:ffv'(ri)...:ffy(r20C'i(t)C„,(t)d^(tKt)}|z), (21) 

fvW^E^Xdn...d..„_, (22) 

{i|r{i?„-(Ti) . . . Hv{T^„-r) Cl{t)d{t)}\i). (23) 
4 Illustration of the method at lowest order 

In order to illustrate the method of calculation, let us consider the first non- 
trivial term {n = 1) in the expansion (21) of . This approximation, 
quadratic in Hy, is formally equivalent to a direct generalization of ND's 
solution of the x-ray problem to several interacting quantum channels a. For 
conceptual details we refer to ND's article [Q. 

The central assumption in this approach consist in considering the localized 
state (j)d to be structureless. The scattering of conduction electrons described 
by H]v therefore depends on the occupancy of the (pd state only, but on no 
other internal structure of (pd- This allows us to consider the scattering po- 
tential Wl^§ as a time-dependent external potential acting on the conduction 
electrons during the times of occupancy of (pd- 

At the order n = 1, for fixed times ti and T2 with, for instance, ti < T2, the 
hybridization term Hy at the time ri describes the tunneling of an electron 
into (pd, whereas the hybridization term at T2 transfers this particle back to 
one of the electrodes. The potential Wj^^ applies therefore during the interval 
{ti,T2) and is zero for all other times. The time t must lie between ri and 
T2 because of the operators d){t)d{t) in the definition of [t). The only 
contribution of the localized state to the amplitude is its free propagation 
factor 

exp(^-iEd{T2 - Ti)). 

If we expand the amplitude {i\T{Hv{Ti)Hy{T2)C}^{t)C'^{t)}\i) in the external 
potential, the Feynman diagrams arising from the different contractions fac- 
torize into two groups: the diagrams that are connected to at least one of the 
points Ti, T2 and t, and the closed loops where no vertex coincides with one of 
these points (see Fig. 2 for an example). If we denote by L"" (ri, T2, t) the for- 
mer group of diagrams, and by e'^^'^'^''^^^ the latter, where by the linked cluster 
theorem the complete closed loop sum can be expressed by the exponential 
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(a) Diagram arising from 
the many-body problem. 



(b) Factorized diagram giv- 
ing rise to three one-body 
problems. 



Fig. 2. Example of the factorization of the many-body problem. The plain lines 
represent free conduction electron propagators (evolving under Hq) and the dashed 
line the propagator of the particle in cp^. The full dots are interaction vertices due 
to W"^" and the open dots at ti and T2 due to F". 

of the sum C(ri, T2) of all single closed loops, the amplitude at n = 1 in (21) 
becomes 

{z\T{Hv{n)Hv{T2)Ciit)C^it)}\t) = L""'(ri,r2,t) e^(---^) e-^^(-^--^) (24) 



As we will see below, both quantities, L°"'(ri, r2, t) and C(ri,r2), can be ex- 
pressed by one-body propagators G"" (r, r') describing a particle which is 
created at the time r' in the electrode a' and annihilated at the time r in the 
electrode a, and which propagates under the external potential W during the 
time interval {ti,T2), 

G-'{r,r') = {^\T{C^{r)CUr')}\^). (25) 
This propagator is determined by the Dyson equation 

G""'(r, = Goi^^ 'r')Saa' r ^0 ^"""'(r", (26) 

a" ^ 

where Go (t, t') is the free propagator. The factor (— i) in front of the integral 
in (26) arises from the unusual definition (25) of the Green's functions. The 
Dyson equation will be calculated using the asymptotic expressions for large 
|r — r'l developed by ND. Since for the free propagators, the electrodes a are 
assumed to be independent, ND's development can be carried out for each Gq 
separately, and is given by (see Eq. (38)) 



V[ I +7rtan79„5(r-r) 
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where the V indicates a principal value distribution. The angle tan da contains 
the short-time corrections to the asymptotic form in l/(r — r') and depends on 
the form of the conduction band. For a completely symmetric band around the 
Fermi surface it vanishes. If we introduce the propagators (27) into the Dyson 
equation (26) and integrate over the delta functions, we find the following 
system of singular integral equations 



i(r)G(r, t') - - r dr" ^(^7^) B{t") G{t'\ t') = /(r, r'). 



(28) 



where the coefficients A and B are 2x2 matrix functions of the form 



A{t) =e^^^A = e^^^(U + 7rtan7?^), 



(^TTfif), 



(29) 
(30) 



for the 2x2 matrices n = {fJ^a^aa'}, tani? = {tsai'&aSaa'}, and g = i/W — 
{i^aVF""^'}, and where 1 is the 2x2 unit matrix. The inhomogeneous term 
/(r, r') is given by 



(31) 



The system of singular integral equations (28) turns out to be exactly solvable. 
We provide the solution in the Appendix A. Let us define the constant matrices 



S ^ A + B = 1 + 7r(tan^? + il)g, 
D ^ A - B = 1 + 7r(tan?? - il)g. 



As shown in Appendix B, the matrix 

e = SD-^ = ~ 7r(tani? + il)g') 



1 -7r(tani?-il)5r 



(32) 
(33) 



(34) 



is the inverse of the (unitary) scattering matrix due to VF""' evaluated on the 
chemical potentials Thus, if we write 



e 



(35) 



the solution of the Dyson equation (26), for r, r' e (^i, T2) is given by Eq. (A. 34), 
and reads 



G{t,t') = i-i) S~\T)e-'^^^^ 







A 




A 


A(t) 






r'-n 


" g+iA(r) 




r-Ti 




r' - T2 





r — T' 



V[ + 7re(T') S{t - t') M{t')u, (36) 
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n t T2 n t 

(a) (b) 

Fig. 3. The two possible contractions contributing to L(Ti,r2,t). The bold lines 
represent the propagator G, given in (36). 

for the matrices 

Q^^'^=U(SD-A)B-' ifr'G(ri,r2), ^^^^ 
Itan-i? otherwise, 

{bD-^B~^ if r' G (ri,r2), , , 

M(r') = <^ ^ ' ^' (38) 

I 1 otherwise, 

S{r) = [^ if-e(n,r.), ^3g^ 
I 1 otherwise, 

A(r) = /^ if-e(n,r2), ^^^^ 
1 otherwise. 

In the limit where the non-diagonal entries of W^"' vanish and the electrodes 
are decoupled, we recover for G"" (Ji't') ND's result ([|ll, Eq. (51)). 

Let us now calculate the factor L"°'(ri, r2, t). Since we know the propagator 
G{t, t') in the presence of W, we will use Wick's theorem to contract the four 
conduction electron operators at the points Ti,T2 and t. Every contraction is 
associated to a propagator G{t,t'). As we see in Fig. 3, there are only two 
possible diagrams. The first diagram. Fig. 3(a), contributes 



-ifY. G'°^"(ri,t)G°'"^(t,r2) 



OL\OL2 



where a contraction of the operators C\t'\C(t^ equals —G{t,t') and where 
the third minus sign arises from the commutation of C'^(t) and Ca'{t) to 
separate the contractions. The second diagram. Fig. 3(b) gives, accordingly, 

(_1)2 J2 G"^"Hri,r2)G'°'°(t,t). 

aia2 

For convenience, we introduce the auxiliary times T3 = r4 = t, and the indices 
as = a' and = a. With these notations we find that 

L""'(ri,r2,t) = Yl detG'"^-^--(r2,^i,r2fe) (41) 

«i«2 



9 



where the determinant is taken over the indices j,k = 1,2. 

If we introduce the exphcit form (36) of the propagator into the latter equation, 
we are lead to divergences in both propagators G{t2, ti) and G{t, t). According 
to ND, these divergences are due to the asymptotic approximation of the 
propagators and must be cured by a (pure imaginary) cutoff i^Q of the order 
of the bandwidth, such that 



Gin,r2) 



1 



n - T2 



-ifJ,{Tl-T2) 



{t^o\r2-n\) - BD-'B-'u, (42) 



and, if we interpret G(t, t) as G{t — i^Q^,t), 

G{t,t) = i-i) S-^BD-^B-^u (i^o)- 



(43) 



Notice that in this form, G{t, t) does not depend explicitly on the boundaries 
Ti and T2. This might be an acceptable approximation for the Green's function 
in the determinant of L"" , but it is not sufficient if we consider the closed 
loop sum. 

The sum over all simple closed loops can be calculated by use of the equation 
g, Eq. (22)) 



A 



dC{n,T2] 



dX 



E 



A=l 



T2 



Pdr Ti [WG{t,t)] 

■J TI 



(44) 



where in C(ri, T2) all potentials W have been multiplied by the parameter A. 

Again G(r, r) is singular, but if we develop G{r, r') in powers of At = t — t', 
we can separate out the singular term which diverges like 1 / At. Namely, since 



X(t) e-'"^^ X-Ht') 



1 


'T-Ti 


# A 
+ — 


"r - Ti" 




r - T2. 


71 


T - T2. 



r -T2 
r ~Ti 



-i/iAt 



r -Ti 
[r - r2) 



T -T2, 



At + C(Ar)2 (45) 



we have in the limit At 



G{t,t') 



Ar^O 



(Sing.) - 



1 1 

+ 



T2- T T - Ti 



S-\t)-M{t)u (46) 

TT 



where the singular part (Sing.) is due only to the approximation on the prop- 
agators. The important point is that it does not depend on ti and T2. As 
shown in ND, its contribution to the exponential is proportional to the first 
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correction in energy (for a static potential W extending in time from — oo to 
+00) which will be denoted by AEi. 



The integration of r over the interval (ti, T2) leads to a logarithmic divergence 
at the boundaries and requires again the introduction of the cutoff so that 
(considering only the regular parts of G(t, t)) 



A 



9A 



A=l 



TT 



log i\T2 - Til^o 



2^ [-1 I.- 
= log i\T2 - Ti 4o 



Tr 
Tr 



o log i\T2 - Ti\Co 



WS-^ABD-^B-^u 

ABD-^S-^ 
A' 



(47) 



The last equality can be shown as follows: If M(A) is an invertible matrix 
depending on the parameter A, we have 



_d_ 
dA 



M-\X) = -M-\\) (^M(A)j M-^(A), 



(48) 



as can easily be verified by deriving the matrix product M(A)M~^(A) = 1 
with respect to A. Therefore, if we replace the potential W by \W in the 
definition of the matrix & in Eq. (34), and derive ©"^ with respect to A by 
using (48), we find after some algebra 

A6^ = -2BD-^S-\ 
dA 

Furthermore, since A is defined by 6 = e^*^, we have 



(-2)IV 



ABD-^S 



ic-l 



ATr 



dA 



= ATr 



Ae^^^ — e-2^^ 
dA 



where we have used the invariance of the trace under cyclic permutations. 
With help of the latter expression we can integrate the simple closed loop 
sum. This yields, reintroducing the contribution in AE'i, 



dC 1 
d^:TV = ^log[i|T2-Ti|eo]TV 

dA TT"' 



Jo 

and the closed loops contribute 



A^ 



- iAEi (t2 - Ti) (49) 



:,C{ti,T2) 



i\r2 -Til^c 



-iAEl{T2—Tl) 



(50) 
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Putting together Eqs. (41) and (50), the contribution to F^' (t) in the series 
(21) takes the form 



(2)' 
\2 



—I 



to<Tl<t<T2<tl 

I dndrs J^iy^'TV"' detG'"^-i"-(r2,_i,r2fc) 



to<Tl<t<T2<tl 



OL\OL2 



^^0^2 - nl 



where rs = r4 = t, as = a', 04 = a, and j, A; = 1, 2. The factor 1/2! occurring 
in (21) has been ehminated by the restriction of the domain of integration. 



5 Extension to higher orders 



An extension to higher orders in perturbation of ND's ideas was carried out 
aheady a long time ago by Yuval and Anderson who derived in this way 
an expression for the partition function of the single spin 1/2 impurity Kondo 
problem. In this section we show how their method can be generalized to the 
present problem. Consider an amplitude at order n > 2 in the Dyson series 
(21) for Fy^ . If the intermediate times are fixed and chosen such that 
Ti < T2 < ■ ■ ■ < r2n, an electron tunnels into the localized state at odd times 
r2j_i, and leaves it at even times r2j (j = 1, . . . , n). As before, these times have 
to be chosen such that the state 0^ is occupied at the time t. The many-body 
problem can be factorized into one-body problems as in the previous section, 
and the amplitude in the expansion (21) can be written in a form similar to 
Eq. (24), 

(i|T{^y(ri) . . . ifv(r2„,) Ci{t) C^,{t)\\i) = 

L""'(ri, . . . ,r2„,t) e^("i'-'^2") g-iB.E;=i (-2.-^2,-1)^ ^52) 



where L""'(ri, . . . , r2„, t) contains all diagrams that are connected to at least 
one of the times ri, . . . , T2n or t, and e'^('^i'---''^2n) jg ^^j^g closed loop sum. Both 
quantities are evaluated in the external potential W{t) which equals W for r 
in the intervals {T2j-i,T2j), j = 1, . . . ,n, and is zero for all other times. The 
one-particle propagator in this potential obeys the Dyson equation 

G(r, r') = Go{t - t') - t T dr"Go(r - r")W^(r")G(r", r'). (53) 

Jto 
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The solution of the latter equation is given in Eq. (A. 40) and reads 



G{t,t') = i-i) S~\T)e-'^^^^ 



-ilJ,{T — T') 



r - 



^+iA(T) 



M{t')v. (54) 



for the same definitions of the matrices as in Eqs. (37)-(39) in which the 
conditions r e (ti, T2) are replaced by r e (''2j-i, T2j), j — 1, . . . ,n. 



In order to determine the factor Z/""'(ri, . . . ,T2„,t) we take advantage from 
the fact that the tj are time ordered, ri < ■ ■ ■ < r2n- A full contraction of the 
operators at the times ti, . . . ,T2n,T2n+i = t,T2n+2 = t can be regarded as a 
permutation tt of n + 1 elements, where the creator C'^(t2j) is contracted with 
the annihilator C(r27r(j)-i)- By Wick's theorem, the whole amplitude must 
be multiphed by the signature of the permutation that brings all contracted 
operators together. Since the operators at the times Ti, . . . ,T2n are explicitly 
time ordered and since the operators at t occur as a pair of creation and an- 
nihilation operators (for which their relative position with respect to the Tj 
is of no importance) this separation is performed by 7r~^. Every contraction 
is expressed by the factor —G°"^^'-^'i-'^'°'^^{T2Tr(j)-i,T2j) as was seen in the previ- 
ous section. Therefore, the sum over all contractions equals the sum over all 
permutation tt and is given by (r2n+i = T2n+2 = t, a2n+i = tt2n+2 = a) 

L°'°' (ti, . . .,T2n,t) = 

n 

{-If Yl (l^"0*'l^"'---'^"'"Esign(7r-')n<^"'^"^-""''(^2.0)-i,T2,) = 

a\...a2n TT j = l 

(-1)" (V^"^)*^"^..V^"^"det,+iG'"^-i'"-(T2,_i,r2fc) (55) 

ai...a2n 



where we used Leibniz' definition of the determinant. The subscript in the 
determinant indicates that it must be taken over the indices j, /c = 1, . . . , n + 
Again the Green's functions G{T2j^i,T2k) and G(t,t) are singular and must be 
cured by the cutoff ^0 as in Eqs. (42) and (43). The closed loops are determined 
by the equation 



A 



dG{Ti, . . .,T2n) 

dX 



= i r dTTr[W{T)G{T,T) 

1 Jto 



(56) 



A similar kind of investigation of the singular and nonsingular behavior of 
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G(t, t) as in Eq. (46) leads to the simple closed loop sum 



n TrA^ 

C(ti, . . . , T2n) = -iA^l E(^2,- - T2,--l) + ^ E (-1)^^' ^^(ieolT,- - T,|) 

i<j 

(57) 

In order to obtain the full expression at order 2n for F^' (t) , we must integrate 
over the times ri , . . . , T2n with the constraint of t lying on a time interval in 
which the localized state is occupied. This can be taken into account if we 
multiply the integrand by the expression Emt=\ (^a2n,.»' (^«2m+i,a '^(t - T2m)5{t - 
T2m+i) and integrate over the 2n + 2 time variables t\ < ■ ■ ■ < T2n+2 and 
sum over the indices ai, . . . , a2n+2- The time order of the integration variables 
drops the factor l/(2n)! in the Dyson series, and the full amplitude reads 



1 



(2n) 



n+l 



= [-1 



dri...dT2n+2 E {V"^)*V"K . .V''^^+^ 



to<Tl<---<T2n+2<il 



«1---Q!2n 



E ^a2m,«' ^oc2m+i,a ^{t-T2m) ^{t-r2m+l) det „+iG"2^-i"2fc (r2j._^, rafe) 



m=l 



exp 



2n 

5^E(-ir^^'K^eo|r,-r,|) 



■ (58) 



6 Discussion and limiting cases 



If wc introduce the expressions (51) and (58) into the Dyson series (21) for 
Fy(^ , we obtain 

oo „ 

(l^")*y"'F-"'(t) = E(-0''^ J dri...dT2„ E 

n p 

E ^«2m-i,a Sa2ra,a' ^(T2m-1 " ^) ^(^ " T2' 

2n 



ai...a2„ 



m=l 



det (^2,-1, T2ik) exp 



^^E(-ir^^l<^eo(r,-r, 

^-i{Ed+^Er) X;j=i(^2j -r2j-i) j-gg^ 



A similar derivation holds for Fy{t). The main difference lies in the fact that 

one single creation operator is fixed at the time t. Its relative position with 
respect to the times of integration is important to determine the signature of 
the permutation separating the contractions. We express this property by mul- 
tiplying the integrand by Zlm=i ^a,a2m.^{'t ~ T2m)- The integration is performed 
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over Ti < ■ ■ ■ < T2n, and the sum over ai, . . . , a2n- This yields 



n=l 



i2n-l 



J2 ^c,,a2m^it -T2r, 



m=l 



dri . . . dT2n 

to<Tl<-<T2n<tl 

(V^'YV' . . . 1/"^" det„G'"2^-i"2'^(r2j_i, 

2n 

t,-a2 1 

exp 



5^^(-ir+Mn(.eo(r,-r,); 



The latter two equations represent the central result of this article. Since they 
are by far too complex for a direct interpretation, we discuss here two limiting 
cases. The first one consists in the scalar problem with only one channel a. 
In the second one we suppress the off-diagonal terms, a 7^ a', in the matrix 



If we assume that the number of quantum channels a is equal to 1, the 
present model can be interpreted as an extension to all orders of perturbation 
of ND's x-ray absorption model [Q, or as the spin 1/2 Kondo problem inves- 
tigated by Yuval and Anderson where the occupancy of the localized state 
(t)d corresponds to the directions of the impurity spin. The propagator at the 
n-th order of perturbation simplifies to 



G(r, t') 



-IV 



n 



T - T2j T - T2J-1 



T - T2j-i r 



V 



T — T 



vrtan 6{ 



T — T 



, (61) 



for r, r' lying on an interval where W ^ 0. The constants (3 and tan-;?' are 
given by (see |l[], Eq. (42)) 



/3 = l- 2TTgtan^ + vr^^ V cos^ 
tan-i?' = tan-i? — ng/ cos^-i?. 



and S is the phase shift. Since all G{T2j-i, T2k) commute and the exponent 5/ti 
is the same for all factors in the determinant, the latter can be expressed in 
terms of Cauchy determinants, which are defined as (see 0, Chap. 7, 3.) 



(Cauchy)^ = det ; 



T2v-1 — T~2k 



n ij2j-l- T2k-l){r2k - T2j) 
^<j<k<n 

n (^2j-l - r2k) 
j,k=l,...,n 



(62) 
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If Ti < • ■ ■ < T2n, we have 

(2n 
E In |r,-T 

Kj 

The determinant in L(ri, . . . ,T2n,t) reads (t2„+i = = 



(63) 



detn+lG{T2j-i,T2k) 



-IV 



n+1 



/5 



3+*Ml]J=l(^2j-T2j-l) 



n+1 

n 



n 

.fc=i 



T"2j-1 - T2k T2j 



T2k-1 



T2j~l — T2fe-1 T2j — T2k 



det 



n+1 



T2j-l - r2fc_ 



where all singularities must be replaced by the cutoff i^o- The last factor is a 
Cauchy determinant in which all factors depending on t cancel out, except the 
singular contribution l/(t — t) for which we introduce a cutoff, too. Similarly, 
the double product in front of the determinant does not depend on t, and is 
nothing else than [i^Q /(Cauchy) Hence 



det n+lG{T2j-l,T2k) 



-IV 



n+1 



^l+2n g+iI]jLi('^2j— r2j_i) 



exp 



- Ti 



i<j 



where we have used the identity J2i<i<j<2n{~^y^'' = If "we define flm = 
< i < j < 2n, ^ {2m — 1,2m)}, we obtain for Fw(t) a result 
similar to that of Yuval and Anderson [01 , 



Fwit) 



11°° 

-— — y 



-IV 



2n 



dri . . . dr2„ 



iO<Tl<-<T2„<tl 

n P n 

E '^('^2m-l - t) 5{t - 1"2m) exp -i{Ed + /^Ei - /i) E(^2j " T-2j-l) 



m=l 



+ (1-2 



<5 , 52 

~r 2 



(64) 



The function Fw{t) does no longer express a current but a joint probability 
amplitude of the presence of a particle in the localized state (pa as well as of a 
particle in the conduction band described by the operator C{t). If we pass to 
imaginary times, tj —i'Tji the latter equation becomes formally equivalent 
to a grand canonical correlation function of a classical one- dimensional gas 
with two different types of particles: a type A situated at the (space) positions 
T2j-i, and a type B at the positions T2j. All particles are subject to external 
potentials, —{Ea + Ai — /i)r for the type A and -\-{Ed + Ai — /i)r for the 
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type B, as well as to logarithmic two-body interactions which are repulsive 
between particles of the same type and attractive otherwise. There is no kinetic 
energy. The cutoff enters as the (hard-core) volume of the particles. In this 
context, the positions T2m-i and T2m that are fixed to the value t by the delta 
functions must be regarded as having the distance between each other. The 
correlation function Fw{t) can therefore be interpreted as the (unnormalized) 
average density of a couple AB at the point t (with A to the left of B). 

We can derive an analogous result for the function Fy{t) which represents now 
the transition rate of an electron into the localized state. In the same way as 
above, we find 



n-1 \ y / to<Ti<-<T2„<ti 
n r n 

^2m) exp -i{Ed + AEi - /i) Xl(^2i - r2j-i] 



'2n 



m=l j=l 

+ (1-2^ + 5) E i-iy^' \n{zUr, - n))]. (65) 

Analogously, this expresses the average number of particles of the type B at 
the position t. 

A investigation of the partition function of the one-dimensional gas was per- 
formed by Anderson, Yuval and Hamann by renormalization group meth- 
ods, a route we will not pursue here. For the general case with the number 
of channels > 2, a similar mapping onto a classical gas is not possible in 
the same way because the determinants of propagators do not form Cauchy 
determinants. Within the scope of the present work a similar, but different 
mapping has not been found. 

Let us now consider the case in which the left and right electrodes are not cou- 
pled via Hw- If we assume further that each electrode contains some number 
of different non-interacting quantum channels, we recover a simplified version 
of the system considered by Matveev and Larkin 0. The latter authors com- 
puted the phase shifts for different geometries of the electrodes and determined 
the current through the resonant level in the barrier in first nontrivial order 
in perturbation, corresponding to [-^\/(^)](2) in our notation. If there is only 
one quantum channel in each electrode, we have 



with a similar definition of as for (3 above. We see that the influence of the 
left electrode enters through the closed loop spreading only. In the limit to ~^ 
—oo, this expression becomes time-independent (it is a Gamma function). 
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whereas scaling arguments reveal the resonant behavior of the current 

[Ed + AEi ^ f^y 



J7.R 



n *B 
• ^~ 2 ? 



(2) 



Following 0] the sign of the exponent depends on whether the impurity state 
is localized near the left or near the right electrode. 



7 Conclusions 



In this paper we studied the tunneling of interacting spinless electrons through 
an insulating barrier. It was assumed that the interaction was restricted to 
scattering of conduction electrons with a single trapped particle in a struc- 
tureless localized state in the barrier. Our central result is given in Eqs. (59) 
and (60) in which we express the tunneling current as an infinite perturbation 
series with respect to the hybridization term of the Hamiltonian. An extension 
of the method of Nozieres and De Dominicis allowed an exact computation 
of each term in this series, including the effects of transient interaction be- 
tween the conduction electrons and the localized state. The renormalization 
of the tunneling amplitude due to the additional charge in the barrier was 
taken into account. 

The calculations were based on two essential assumptions: a separable po- 
tential, and a completely structureless localized state. Since the uncoupled 
systems were assumed to be Fermi liquids, the model does not apply to one- 
dimensional systems. With these limitations, extensions of the model may be 
considered, such as to include electron-phonon interactions, many traps, or 
disorder in the barrier as well as in the electrodes. Because of the nanometric 
size of the system, disorder in the electrodes situates us in the regime of the 
universal conductance fluctuations (UCF) which have been intensively inves- 
tigated in the last two decades [0. It is beheved that further developments 
of the present work will yield a microscopic description of the stochastic time 
behavior of these fluctuations. 



A Solution of the singular integral equation 



In this appendix we solve the system of coupled singular integral equations 
obtained from the Dyson equations (28) and (53) for the one-body matrix 
propagator (r, r') which we denote in this section as yj"" (r, r'). Generally, 
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these systems of integral equations can be written in the form 



A(rMry) - 1 / £(l!MlVO dr» = /(r.r'). (A^l) 

Here A, B and / are given N x N matrix functions, is the unknown prop- 
agator with the entries , and L is a contour in the complex plane with 
the end points t and t' which may take finite values or lie at infinity (we do 
not consider the case where L is a closed contour, i.e. t = t'). The positive 
direction on L is chosen to run from t to t'. 

The technique to solve such integral equations was developed by Muskhel- 
ishvili [^] and Vekua ^ and consists normally in a transformation to a system 
of ordinary coupled Fredholm equations. However, the particular form of the 
coefficients A and B as given in Eqs. (29) and (30) and their generalization 
to higher orders in perturbation allows us to develop a compact, closed so- 
lution by a method similar to that used by Yuval and Anderson for the 
scalar problem. This method is in fact nothing else than Vekua's technique to 
regularize discontinuities in the coefficients A and B, and we formulate it for 
arbitrary smooth contours L in the complex plane. In Sec. A. 3 this solution 
will be carried out in detail. Furthermore, we will show its limitation in the 
solution of more general problems. 



A.l Transformation to a Hilbert problem 



The first step in the solution of Eq. (A.l) consists in a transformation of the 
integral equation to a Hilbert boundary value problem. We follow the lines 
of Muskhelishvili and Vekua (see e.g. §16), which we expose briefly in the 
sequel. 

From now on it is assumed that the functions A{t) and B{t) are N x N 
matrices, which are - for the theory to be valid - required to be such that 

S{t) ^ A{t) + B{t), D{t) ^ A{t) - B{r), 

are regular, i.e. det 5'(r) ^ and detD(r) ^ for all t E L. Furthermore we 
assume that all entries of A{t) and B{t) satisfy a Holder condition on the line 
L except at some points ri, . . . , G L where they have flnite discontinuities, 
but where they still are assumed to be such that S{t) and D{t) remain regular 
when considering the left and right limits r ^ ± {k = 1, . . . ,d) on L. 

The Hilbert problem is formulated by means of the function \l/(-z), which is 
holomorphic in the whole plane, except on the line L where it is continuous 
from the left and the right side with respect to the positive direction of L 
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(except possibly at the points r^) 



TTZ JL T" — Z 

In the terminology of Muskhelishvili and Vekua, such a function is called a 
sectionally holomorphic function. 

The connection between the functions \& and cp is given by the Plemelj formulae 
(see e.g. §, §17or i,§3) 



v[/+(r)-vl/-(r)] =fi(r)yp(r,r'), 

1 r B{T")ip{T",r') 



^+(r) + ^"(r) 



m JL 



t" -T 



dr" 



(A.2) 



for T & L, not a point of discontinuity r^. The functions \l/"'"(r) and ^I/~(t) 
denote respectively the limiting values of "^{z) when z approaches the point 
T & L from the left or the right side of L with respect to the positive direction 
on L. By replacing the integral in equation (A.l) by the second Plemelj formula 
(A.2) we obtain 



A(r)y,(r) = l[vl>+(r) + vl/-(r) 



fir). 



Adding and subtracting the first Plemelj formula (A.2) from this last expres- 
sion leads to 



v{r) = [S{T)]-'^^T) + [S{Tr'f{T), 
v{r) = [D{r)]-'^-{r) + [Dir)]-'fir). 



(A.3) 



From these two expressions we obtain the following Hilbert boundary value 
problem, equivalent to the integral equation (A.l), 



^+(r) 



where we have defined 



6(r) = Sir) [Dir)]-' = [Air) + i?(r)] [Air) - fi(r)]-\ 
bir) = (&ir)-l)fir,r'). 



(A.4) 



(A.5) 
(A.6) 



The matrix (5 is called the coefficient of the Hilbert problem. Once we have 
found the solution of the Hilbert problem (A.4) the function ifir, r') can be 
recovered by any of the equations (A.3). 

As will be shown below, this matrix (3 is equal to the inverse scattering matrix 
at the Fermi levels, and is therefore unitary. 
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A. 2 Formal solution of the integral equation 



The Hilbert problem (A. 4) is solved by means of a particular solution (the 
so-called fundamental solution) of the homogeneous Hilbert problem obtained 
by putting 6(r) = in (A. 4). This fundamental solution will be computed in 
Sec. A. 3 and is given by a x matrix function X[z), where z is a point 
in the plane, not on L. The limits of X{z) on L from the left and from the 
right with respect to the positive direction on L will be denoted by X~^{t) 
and X^(r), respectively, for r G L\{ti, . . . , r^}. 

Since X{z) is a solution of the homogeneous Hilbert problem corresponding 
to (A. 4), we have 

G{T)=X^r)[X-{r)]-\ 

By Vekua 0, §6, the general solution of the non homogeneous problem (A. 4), 
having at most polynomial divergence at infinity, is then given by the expres- 
sion 

X{z) f [X+{t")]-H{t"'' 



^(^) = ^ / ^ hr dr" + X{z) p{z) 

with p{z) a matrix of size N x N with arbitrary polynomials as entries which 
must be determined by the boundary conditions. In the physical applications 
to follow, it is sufficient to notice that the function ip{T, r') is a Green's function 
which must vanish as its arguments go to infinity. Since the function X{z) 
remains bounded at infinity (as we will see later) these polynomials must 
identically vanish, which allows us to drop the term X{z)p{z) in the sequel. 

The boundary values \E'~''(r) and ^E'^(t) follow again from the Plemelj formulae 
and are given by 



The solution of the original integral equation (A.l) can then be computed by 
any of the equations (A. 3) and takes the form 

1 r B(t") fir") 

V.(r, r') = A{t) fir) - - Z{r) / ^ ^ dr" (A.7) 

m Jl t' — t 
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Fig. A.l. Contour C for the integration of formula (A. 11). 



where we have defined 



Ar) = \{[S{r)]-^ + [D{T)]-'], (A.8) 

B{T) = \[[X\T)r-[X-{T)]-'Y (A.9) 
Z{t) = [S{t)]-' X^{r) = [D{t)]-' X-{t). (A.IO) 

The fact that the inhomogeneous term f{T,T') is proportional to the free 
propagator and that B{t) is the difference between the hmiting values of the 
inverse fundamental solution allows us to evaluation the remaining integral in 
(A.7), 

Jl t" — T 



Let us write 



where a and b are N x N matrices that may depend on r' but not on r. Such 
/ can be continued analytically to the remaining complex plane, and vanish 
at \z\ — > oo. Considering the term proportional to a only and assuming that 
T 7^ r', the integral (A. 11) can be calculated explicitly by the integration of 
[X(z)]'~^a/ (^(z — r)(z — r')^ over the contour C indicated in Fig. A.l. In fact, 
by the theorem of residues, we have on the one hand 

/ [X{z)]'Udz _^ 
Ic {z — t){z — t') 

Since X{z) is the fundamental solution, we will see in Sec. A. 3 that its entries 
are non- vanishing, analytic functions everywhere in C except on L\{ri, . . . , Td] 
where it is continuous from the left and the right. At the points Ti,...,Tii 
it either has integrable poles or vanishes. Moreover, we will see that X{z) 
remains bounded at infinity, and hence there is no residue at infinity. On the 
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other hand, the integral on the contour C can be spht as follows 

[X{z)]-^adz f [X+{T")]-^adT" f [X-{t")]'' adr" 



c {z-t){z- t') Jl {t" - t) {t" - t') Jl {t" - t) {t" - t') 

— m : a — m : a. 



T — T 



T — T 



Notice that the singularities at the points of discontinuity (if there are any) 
are of order z~'^ for certain < £ < 1 and do not contribute. When combining 
the last two equations, we obtain the integral, 

- {[XHt')]-' + [X-(r')]-') « 

+ ^(i^"Mi-'-ix-(T'r ')(-}. 

where we have now added the contribution arising from the integration over the 
delta function. Since we assumed explicitly that r ^ r' the last expression has 
to be considered as being a Cauchy principal value distribution. The singular 
contribution for r — r' can be computed from the Poincare-Bertrand formula 



V\ 



1 



V 



1 



V\ 



1 



T — T 
2 Xl^'l 



V 



1 



1 



+ 7r'5(r" - t)5{t" - t') 



and is given by 



2 

\{[X^{T)]-'-[X-{T)]-')a5{T-T'). 

These results are introduced into the expression (A. 7) which leads to the 
solution of the integral equation (A.l), 



^{t,t')^Z{t)[Z{t')]-' {A{T')a + ^.C{T')h) ^(:^) 
+ {A{T')h + mC{T')a) 5{t-t') 

where A{t) and Z{t) are defined by (A. 8) and (A. 10), where 



(A.13) 



C{r) = - { S{t) 



D(t] 



(A. 14) 
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A. 3 Fundamental solution of the homogeneous Hilbert problem obtained from 
the integral equation 

In this section we determine the fundamental solution of the homogeneous 
Hilbert problem 

<|.+(r) = 6(r)<l>-(r) (A.15) 

for r G L\{ti, . . . , r^}. The matrix &{t) is defined by Eq. (A. 5) for the differ- 
ent cases imposed by the underlying physical problems. 

A. 3.1 Characterization of the fundamental matrix 

Following Vekua 0, §5, a fundamental solution of (A.15) is a matrix, X{z), 
characterized by the following two properties: 

(a) If X{z) is a fundamental matrix, then its determinant, detX{z), does 
not vanish anywhere in the finite part of the complex plane. 

(b) Let ((z) be the vector formed from the column /?(/? = 1,...,A^) of 
X{z), Cai^) = Xap{z) for a = 1, . . . , A^, and (— k/j) its degree at infinity 
(i.e. Ca{^) — '^'^) -2 ^ oo). Then the determinant dety(2;) of 

/3 

the matrix Yap[z) = z'^^C^^z) has a finite non-zero value at infinity. 

Any system of solutions (represented by such a matrix X) satisfying these two 
properties will be called fundamental, and it can be shown (0, §5) that the 
set of degrees at infinity . . . , k^v} is the same for all fundamental systems. 

A. 3. 2 Regularization of the discontinuities 

Consider the homogeneous Hilbert problem (A.15). We assume first that the 
matrix &{t) is a constant, unitary matrix (3 on the whole line L, G{t) = 
6,Vr G L. 

We choose an auxiliary line L' in the complex plane such that C = L U L' 
forms a simple, closed, smooth contour with the positive direction running 
from t to t' on L and from t' to t on L'. The domain bounded by C, denoted 
by V~^, is assumed to lie to the left of C with respect to the positive direction. 
The unbounded region outside the contour is denoted by V~ . As long as these 
requirements are fulfilled, the choice of L' is arbitrary and has no influence on 
the result. The Hilbert problem (A.15) can be extended to the contour C by 
putting ©(t) = 1, the N x N unit matrix, on L'. With this choice ©(t) has 
discontinuities at the points t and f. 



24 



The method presented below consist in a simpUfied version of Vckua's tech- 
nique to suppress the discontinuities of ©(t). Generally, this transforms the 
coefficient 6(t) into a smooth function only and does not give a solution to 
(A. 15). However, in the present case of a constant, unitary coefficient it is 
already sufficient to solve the problem. 



-1 



6(^ + 0) e(i-o) = e 



-1 



e{t' + 0) e{t' - 0) = e. 



(A. 16) 



The unitarity of & assures that we can find a unitary matrix Bq that diago- 
nalizes simultaneously 7 and 7' such that 



X^Bo'-fBo, 
A' = So-S'So, 



(A.17) 



are diagonal and inverse to each other. Further we may choose the real numbers 
5a, with — TT < 5a < +7r, a — 1, . . . , N, such that 



A« 
A' 



-2i5c, 



(A.18) 



The condition — tt < 5a < +7r is a mathematical requirement that ensures the 
convergence of the expressions, but it is shown in Appendix C that the same 
restrictions result from purely physical arguments. Even more, the physical 
constraints fix uniquely the angles 5a to be obtained from the phase shifts 
from the quotient of the scattering matrices immediately before and after a 
transition (see Eq. (A. 16)). 

We choose an arbitrary point Zq in T>~^ and define the diagonal matrix functions 

^0, and ^ as 



We define the matrices 

7 = 
7' = 



(Co(^))^ = - Zo)'- , 



(A.19) 
(A.20) 

(A.21) 



for a = 1, . . . , iV. Let £ be a not self- intersecting smooth line that starts at 
Zq, intersects C once in t and goes then to infinity. The whole fine i is chosen 
to be the branch line of ^o(^); its part starting at t the branch line of ^i{z), 
and its segment {zQ,t) the branch line of ^{z). Similarly we construct the 
functions ^q{z),^[{z) and Ci^) which are given by Eqs. (A.19)-(A.21) with 
the substitutions t' for t, and with a branch line i' that starts at Zq, crosses C 
at t' and goes then to infinity. 



25 



In this way, the functions ^{z) and ^ire holomorphic in V~ , whereas 

^i{z) and Cii^) ^re holomorphic in V^. All four functions can be continued 
analytically to the contour C with the exception of their respective branch 
points t or t'. 



Now introduce the matrices 



Ao = &{t + 0) Bo [Ut + 0)] ' = &it - 0) Bo [Ut - 0) 
A', ^e{t' + 0) Bo [Ut' + 0)] = 6(i' - 0) Bo [Ut' - 0) 



(A.22) 



and define the functions 



Q(z)^AoUz), 
Q'iz)^[Q{t')]''AU[{z), 

for z e 1?+, and 

P{z) = Boaz), 
P'{z) = [P{t')]''Boeiz), 

for z e V". These functions allow us to substitute the sectionally holomor- 
phic function ^'(-z) in the homogeneous Hilbert problem (A. 15) by the new 
functions 





Q'iz) 


-1 




< 




-1 






'P\z) 




'P{z) 



^{z) for z e 
^{z) for z eV~, 



(A.23) 



which leads to the new Hilbert problem 

with the new coefficient 

6(r)=[g'(r)]"' [Q{T)]-'e{T)P{T)P\T). 



(A.24) 



(A.25) 



An investigation of the phases of the six functions S,...{'t) and ^' (r) shows that 
(3(r) equals the unit matrix on the whole contour C, including the points t 
and t'. The Hilbert problem (A. 15) has therefore been reduced to the trivial 
problem ^+(t) = ^~(t), which consist in seeking a function which is holo- 
morphic in the whole plane C. This condition, together with the requirements 
of Sec. A. 3.1, shows that the fundamental solution of (A.24) must be a con- 
stant matrix, which we choose to be the unit matrix, X{z) = 1. Inverting 
the relations (A.23) leads to the fundamental solution X(z) of the original 
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problem (A. 15), 



^ ' ^ P{z)P'{z) iorzeV-, ^ ' 

and it can readily be verified that this function also satisfies the requirements 
of Sec. A.3.1. 



A. 3. 3 Extension to a coefficient fluctuating between two values 

Let us now extend the above method to the case where 6(t) is of the form 



6(r) 



6 for re (Tj,Tj), j^l...n 
1 otherwise 



where ri, r{, r2, . . . , r„, are subsequent points (with respect to the positive 
direction) on L, and where 6 is a unitary matrix. 

We now apply the same regularization method as above. Let us introduce 

k k fc k k k 

the functions ^o(-^))^i(^) ^{z), as well as ^' o{z) , ^' i{z) and ^^i k — 

1, . . . ,n, which are defined by Eqs. (A.19)-(A.21) with the substitutions Tk for 
t and for t', respectively, and the appropriate choice of their branch lines. 
Consider now the points ri and t[ and define 



P'iz) 



Q{z) = AoUz), 
Q'{z)^\Q{T[)r'A',l{z), 



where 



io = e(Ti±o)So[eo(Ti±o))] \ 

A',^e{T[±Q)Bo[i',{T[±0))]~\ 
The new coefficient of the Hilbert problem is given by 

e«(r) = [Q'(t)]~' [Q{r)]-' e(r) P(r)P'(r). 

As we have seen above, in the interval {ti,t[) this new coefficient is equal to 
1. But since &{t) takes only the two values 1 and 6, which are both unitary 
matrices, and since Bq diagonalizes 6(t) on the whole line L, the module of 
is equal to 1 everywhere on L. However, the phases vary between the 
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different segments on tfie line, sudi tliat 

,2iA £qj, ^ ^ (ti,t() and wliere 6(r) = 6 



e«(r) 



otlierwise 



for the diagonal real matrix A = diag {5i, . . . ,5f^). 
If the original Hilbert problem is written in the form 

^+(r) = 6(t) *-(t) 
we have now found an equivalent Hilbert problem expressed by 



C+(t) 



vl>-(r) 



where 



C+(r) = g(r) g'(r), C-(r) = P(r) P'(r) 



This is a system of scalar (decoupled) Hilbert problems with discontinuities 
at the points T2, Tg, . . . , t„, which can readily be solved. We can use the 

k 

functions ^._o,i(''') eliminate the jumps of ©^^^ between the segments. What 
remains is again a trivial Hilbert problem, ^+ = for which a fundamental 
solution is given by X = 1. In fact, we introduce into "^{z) for each point of 
discontinuity Tk the factor 



k k 



Ur'k)kK + 0) 

1 -1 



for z e P+, 
for z e V~. 



Notice the factors with rj^ as argument. Their effect is to assure the validity of 
the solution to the Hilbert problem for any choice of the origin of the phases 
of complex numbers (i.e. the axis with respect to which they are measured). 
The final fundamental solution to the Hilbert problem, denoted by X[z), can 
then be written in the form 



X{z) = 



fc k 

Bo Y[Uz)^\{z)B^' 



for z eV+ 

Bo n e(^) e'(^) Urk) Coirk + 0) ^o"' iov z e V 



k=l 

" fe k k 



(A.27) 



k=l 



where we have chosen to take {o(''"fc + 0) i^^o the definition of "if (z) and we 
have multiplied all expressions by Bq^ from the right (which does not change 
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the Hilbert problem, but assures that X{z) tends to the unit matrix 1 as 2; 
goes to infinity). 



A. 3. 4 Limits of the method 



We show briefly why the above method of regularization does not give the 
solution for the case where &{t) is a more general piecewise constant function, 
taking the values 61, 62, ■ ■ ■ on L. 

For instance, assume that the coefficient has the form 

'61 for T e /i = (ti,T2) 
^ \&2 for T e 72 = (t2, T3) 
1 otherwise 

for Ti,T2, Ta some subsequent points on L. 

For the Hilbert problem 

^+{t) = e{T)-^-{T) (A.28) 

the regularization above consists in an ansatz in form of a product, "^{z) — 
■01 (-2) '02 (-2) where ipk{z) is the solution of 



&k i^k (t) ifrelk 



il^j. (r) otherwise 
for k — 1,2. If we put this expression into Eq. (A.28) we must have 

i^t{r)^Pt{r) = &{T)i^^{T)^P^{T). 

This equation is satisfied for t & Ii. But if r e I2, the function ■02 (-2) must 
simultaneously satisfy the equations 



iIj+{t) = 62^2 



In other words, since iPi{t) = ipi (r) = '0i(t) on I2, it is required that 



This means that 62 and '0i(r) must commute, which is equivalent to saying 
(by the construction of ipi) that 62 and &i must commute. Generally, there 
is no reason why they should, and the regularization method sketched above 
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cannot be used to solve the problem. In such cases, Vekua's complete theory 
which uses Fredholm integral equations is required. 

If, however, &i and ©2 (or all matrices in an extension to more than two 
different matrices) commute, they can all be diagonalized simultaneously, and 
the procedure of the previous section applies. The only modification in the 
result (A. 27) is that the matrices A in the definition of the functions ^ have 
additional indices k indicating their origin from [&k]^^ Gk-i- 



A. 4 Solution of the Dyson equations 



We turn now to the solution of the Dyson equations (28) and (53). Consider 
first Eq. (28) which is almost of the form of the singular integral equation 
(A.l) considered in the previous section. The only difference in (28) is the 
dependence of the coefficients A{t) and B{t) on the phase factors e^^^. In 
the transformation to a Hilbert problem this difference is of no significance. 
Therefore, if we follow the lines indicated in Sec. A.l we are lead to the 
equation 



(A.29) 



with 



-1 



6(t) = [A(r) + B(t)\ [A(r) - B{t) 

'A{t) + B{t)\ [A{t) - B{t) 
6(t) e-^''^ 



= e 



and 



6(r) = (6(r)-l)/(r,r'). 



The coefficient &{t) fulfills the requirement of being a constant unitary matrix 
function on the line L = {ti,T2) with the value (from Eq. (34)) 



6 = (U - 7r(tani? + il)g^ 1 - 7r(tani? - il)g 



(A.30) 



By comparison with the result of Appendix B, Eq. (B.8), this function is the 
inverse scattering matrix for electrons at the Fermi surfaces, and we can define 
the hermitian matrix A by 



6 



(A.31) 



The eigenvalues of A are the phase shifts which we choose, as discussed in 
Appendix C, to lie between — tt and +7r. 
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Hence, the homogeneous Hilbert problem 



$+(r) = 6(r) $-(t) 

obtained from (A. 29) by the substitution ^(z) — e^^^^^{z) and by putting 
6(r) = 0, can be solved by the method described above. This leads to the 
fundamental matrix X{z), given by Eq. (A. 27), 



X{z) 



B,Uz)i[{z)B^' 



for z e P+, 



Bo i'{z) Ur',) + 0) i?o"' for z G P". 



Since the line L coincides with the segment [ti, T2] of the real axis, this formula 
can be simphfied. If we choose for the logarithm the usual cut on the negative 
real axis, the function X[z) can be written in the form 



X{z) 



Z-T2 



Bn 



Z-T2 



Bo' 



.Z-TxJ \Z-T2, 

for any not on [ti,T2], with the definition — exp(Alogi). 



(A.32) 



To establish the fundamental solution, X{z\ of the original homogeneous 
problem obtained from (A. 29), it is not sufficient to^substitute X{z) = e*^^ X{z), 
because this function is singular at infinity. For X{z) to be bounded at infin- 
ity, we use the fact that X{z) — > 1 as |z| — > 00, and that a solution to the 
homogeneous Hilbert problem remains a solution when it is multiplied by any 
matrix from the right. Therefore, the matrix function 



X{z) 



X{z) e 



■l/JyZ 



(A.33) 



fulfills all requirements for a fundamental solution to the original homogeneous 
Hilbert problem. 

The inhomogeneous term / is given by Eq. (31), 



f{T, t') = e^^^ Go{t - t') = i-iu) V[ + TT tan 1? 5{t - r') 



T — T 



and is thus of the form of Eq. (A. 12). The solution of the integral equation 
(28) is then provided by Eq. (A. 13), which reads explicitly, for r, r' e M, 



G{t,t') = {-i)S-\T) e-^^(^) 



T-T2 



r -Ti 



Tl 



T2 



^(7377) + ^^(^0 - ^0 M{T')iy, (A.34) 
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where we have used the identities 



U -itaiii? = {S-t)B-\ 
l + itani? = {1-D)B-^, 



and where we have defined 

e(r') = 

M(t') = 
A(r) = 



«(5L>->l)S-i ifr'e (ti,T2), 
tan 1? otherwise, 

BD-^B-^ if r' e (ri,r2), 
1 otherwise, 

S if r G (ri,r2), 
1 otherwise, 

A ifrG(ri,T2), 
otherwise. 



(A.35) 
(A.36) 
(A.37) 
(A.38) 



Now consider the Dyson equation (53) at higher orders in perturbation. The 
new feature is that the potential W is nonzero only within the time intervals 
{T2j-i,T2j), j = l,...,n, on the real axis. Hence only on these intervals the 
matrix & takes the value (A. 30), everywhere else it is equal to the unit matrix. 
By Eq. (A. 32) the fundamental matrix reads 



^(-)= n 



Z - T2j-l 



^0 n 



'2j 



Z - T2,_i 



-^0 ^ 



(A.39) 



and this expression must now replace the in the solution (A. 34), 



G{t,t') = {-i)S-\T) e-^^(^) 



T - T2i 



r - T2,.i 



^-inir-T') g+iA(r') 



7=1 ^ 



T — T 



- ) + nQir') S{t - r') M{t')u, (A.40) 



with corresponding definitions of the matrices Q{t), M{t), S{t) and A(r) as 
in Eqs. (A.35)-(A.38). 



B The matrix & is unitary 



We show that the matrix (3 appearing in the Hilbert problem obtained from 
the Dyson equations is equal to the scattering matrix S associated to the 
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external potential W. This is a generalization of a result of Kohn [|T0| who 
determined the phase shift by similar scattering theory arguments for the 
scalar case. In this appendix Hq denotes the free Hamiltonian and H = Hq+W 
the Hamiltonian with the static potential. 



B.l T operator and S matrix 

Let G{z) and Go{z) be the resolvents 



1 

z-H , 



z-Hn 



G 
Go 

for a point z in the complex plane, not in the spectrum of H oi Hq. 
The operator T{z) is defined as follows 

T{z) =W + WG{z)W, (B.l) 

and satisfies the two identities 

Go{z)T{z) = G{z)W, 

T{z)Go{z) = WG{z). ^ ' 

Combining the definition (B.l) with (B.2), the Lippmann-Schwinger equation 
for the T operator follows 

T{z) = W + WGo{z)T{z). (B.3) 

In a plane wave basis of electrons with momentum p and additional quantum 
number a, the scattering matrix S is connected to the T operator by the 
relation (see [|ll|]) 

{p'a'\S\pa) = 6^^\p — p') 5aa' — 2TTi6{Ep/a' — Epa) lim {p'a'\T{Ep + ie)\pa). 

(B.4) 

B.2 Case of a separable potential 



Consider now the physical situation described in Sec. 2. The potential W is 
assumed to be due to a localized scatterer and sufficiently short ranged. Its 
matrix elements between the plane wave states \pa) can then be written in 
the form (see Eq. (5)) 

V V ppl ujp V V ujpi 
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with 

/ dp|ti"P = 1 and u" = 1 

for any vector on the "Fermi surface" of the channel a. With this envelope 
function we can define the orthonormal localized conduction electron state 

= Jdpu^lpa) 
in which the potential decomposes as 

W = Y,\C)W"''' {C'\. (B.5) 

aa' 

Taking the matrix elements between and \C"') of the Lippmann-Schwinger 
equation (B.3) yields 

P 

where T""' = (C"|T|C"'), and 

,/3|2 



G^,{z) = {C^\Go{z)\C^) = Jdp- 



z — Epfj 

In matrix notation for the indices (a, a') this becomes 

T{z)^W + WGo{z)T{z) 

which has the solution 

T{z) = [1-W Go{z)]-^W = W[1- Go{z) W]'^. (B.6) 

An expression for the S matrix (seen as a finite matrix with indices (a, a')) 
can be found from (B.4) together with (B.6). This yields 

{p'a'\S\pa) = 6^^\p - p') 6aa' - 2m d{Ep^ - E^,^,) lim {p'a'\T{Ep^, + ie)\pa) 
= 5^^\p - p') 5oca' - 2ni6{Epa - Ep,^>) 

= b^^ {p - p') S^a' - 2nid{Ep^ - Ep,a') 

u-^: (W[l - Go{Epa + iO)W]-Y" <*■ 
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An integration over the momenta p and p' for the fixed energies E and E' then 
gives 

S'''''{E,E') = J dpdp'6{E' - Ep,^,)6{E - E^^) {p'a'\S\pa) 

= DJaa' - 2m u""' (E) (w[l - Go{E + iO)W]-^ 



u'^iE) V^{E) 



6{E - E') 



= 6{E-E'), 

where Va{E) is the density of states at the total energy E of the system. 

We fix the energy of the incoming electrons at the value /Xq, and consider a 
scattering onto an outgoing state with energy /Zq,/. Since generally /Iq, 7^ /Xq,', 
the energy difference must be contributed by the localized scatterer, such that 
E = fia + Eis^a = fJ'a' + -Eis.o' remains constant, for -Eis,o the initial and E\s^a' 
the final energy of the scatterer. This allows us to calculate 



G^{E + iO) = / dp 



\u. 



a\2 



E - Epa + iO 



for epa = Epa - Eis,a- If Va{e) = z/(e + 
value of the integral, 



dp 



\u. 



a\2 



and if V denotes the principal 



G^{E + iO) = Jde 



\u?\''uJe) 



/ia — e + iO 



a |2 



+ ?de 



Ha - e 



= —Txi — 1^ tan-i^a 
where fa = t'a(Aia)- The angle -da is given by the relation (see also [Q], Eq. (39)) 

1 



tan 'dr 



de 



e - /i„ 

and is identical to the angle appearing in Eq. (27). Hence, 

5"'"(E) = z/„ 6aa' - 27ri Ua' {w[l + 7r(tan^9 + it)iyW]'^ 



or, in matrix notation with 1/ = {I'aSaa'}, tan^} = {tan-f^a 5q,q,/}, and g 
S{E) = (1 - 2ni g[l + 7r{tan^ + il)g]-^)u 



1 + TT^tani!) - il)g) 1 + 7r{tan'd + il)g v (B.7) 
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B.3 Unitarity of the matrix & 



In order to check the unitarity of the matrix given in (B.7), we investigate 
the product S{E)*S{E'): 

= ^ / dpdp'dp"S{E,^ - E) 5{E,,^, - E') S^:*S^^:'. 

a" a" ^ 

The unitarity of S{p, p') imphes that 

E / dp" s^! * s^^: = b^^, {p - p') 

a." ■' 

and hence 

a" 

= y dp dp' - E) 5{Ep,a, - E') 5(3) (p - p') = 5(E - E'). 
In other words 

S{EyS{E') = vH{E-E'). 

Therefore, the matrix S{E)v^^, expressing the scattering matrix at fixed 
energy E per state at the Fermi surfaces, is unitary By comparison with 
Eq. (A. 30) we see that 

e-^ = S(E)u-^ = (1 + TT (tan 1? - il) uW) [l + n (tan 1? + il) uW] (B.8) 

which is the same matrix that was obtained in the Hilbert problem of Sec. 4 
and 5 and yields the physical interpretation of this matrix &. 



C Discussion of the choice of the phase shifts 



In Sec. A. 3. 2 we have seen that the numbers 5a are not uniquely defined, and 
can be chosen to be positive or negative within the limits ±1. It will now be 
shown, that the physical constraints determine uniquely all these numbers. 

For this purpose, let us first focus on the scalar case (A^ = 1). The matrix 
&{t) — e^*''^'^) is then a scalar, piecewise constant function, where S{t) is a 
real number. As has been shown in Sec. B, 6(r) is nothing else than the 
scattering matrix on the Fermi surface. This explains why ©(r) is unitary, 
and that the values 5(r) are the phase shifts which the wave functions acquire 
during the scattering on the local potential W{t). These phase shifts are 
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Go 



1 




i 



& 



d-l 



n T2 T3 

Fig. C.l. A possible form of &{t) 



defined in [— tt, +7r] modulo 2% only, and for two different values of &{t) on 
L the offset from this central interval, 27rm for m an integer, is generally not 
identical. A possible situation is sketched in Fig. C.l. The arrows in this figure 
indicate the increments or decrements of 5(r) at the points of discontinuity. 
Exactly these in-/decrements enter the exponents of the functions ^ defined in 
(A.19)-(A.21). The condition — tt < S < +7r tells us moreover that the different 
phase shifts S{t) can only differ by less than ±27r. They must therefore all have 
the same offset 27rm from the interval [— tt, +7r]. As long as this condition is 
satisfied, the value of m is of irrelevant. For instance we may take m = 
(which is perhaps the most intuitive choice for the regions where S = 1). 

The main condition is therefore seen to be the requirement — tt < Sa < +vr 
but it must be emphasized, that it is so far a purely mathematical constraint 
which ensures the integrability of the equations. However, there is a physical 
argument for this choice which was given in ND's paper: only the perturbative 
solution must be retained. Take Eq. (A. 40) for G{t,t') into consideration. We 
know that the function contains the factors 



for the different 5k — S{Tk + 0) — S{Tk — 0) at the points of discontinuity r^. If 
we let the scattering potentials tend to zero, H^(r) — ^ 0, we must obtain the 
free propagator G^^\t — t') given in (27). If in this limit 5{Tk + 0) — >• 27rm+ 
and 5(Tfe — 0) — > 27rm_ (for m+,m_ some integers), there would still remain 
the factor 



2(n+— n_) 



in G{t,t') which does not appear in G^°\t — r'). Hence the requirement 
m+ = m_ = m at all Th- 
in the matrix case where all matrices &{t) commute (for any other case, we 
do not possess a closed solution, as has been shown in Sec. A. 3.4) all matrices 
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can be diagonalized simultaneously. The system transforms then into a set of 
decoupled scalar problems to which the same arguments as above apply. 



Acknowledgements 

We thank B. Altshuler, J. -P. Ansermet, H. Kunz, P.A. Lee, P. Nozieres and 
T.M. Rice for stimulating discussions. The support of the Swiss National Fonds 
is greatfuUy acknowledged. 



References 

[I] P. Nozieres and C. T. De Dominicis, Phys. Rev. 178 (1969) 1097 
[2] P. W. Anderson, Phys. Rev. Lett. 18 (1967) 1049 

[3] G. Yuval and P. W. Anderson, Phys. Rev. B 1 (1970) 1522-1528 

[4] P.W. Anderson, G. Yuval, and D. R. Hamann, Phys. Rev. B 1 (1970) 

4464 

[5] K. A. Matveev and A. I. Larkin, Phys. Rev. B 46 (1992) 15337 
[6] L. V. Keldysh, Zh. Exsp. Teor. Fiz. 47 (1964) 1018 [Soviet JETP 20 
(1965) 1018] 

[7] G. Polya and G. Szcgo, Aufgaben und Lehrsdtze aus der Analysis, Vol. II. 

Springer- Verlag, 1971 
[8] N. I. Muskhelishvili, Singular Integral Equations. P. Noordhoff Ltd., 

Groningen, NL, 1953 
[9] N. P. Vekua, Systems of Singular Integral Equations. P. Noordhoff Ltd., 

Groningen, NL, 1967 
[10] W. Kohn, Phys. Rev. 84 (1951) 495 

[II] J. R. Taylor, Scattering Theory: The Theory of Non-relativistic Collisions. 
John Wiley & Sons, New York, 1983 

[12] For reviews see: P. A. Lee and T. V. Ramakrishnan, Rev. Mod. Phys. 57 
(1985) 287; A. J. Efros and M. Pollak (eds.). Electron- Electron Interac- 
tions in Disordered Solids, Elsevier, 1985; B. L. Altshuler, P. A. Lee, and 
R. A. Webb (eds.), Mesoscopic Phenomena in Solids, Elsevier, 1991; and 
references therein 



38 



